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ABSTRACT 

o : 
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It seems necessary to suppress, at least partially, the formation of structure on 
subgalactic scales. As an alternative to warm or collisional dark matter, I postulate 
^ | a condensate of massive bosons interacting via a repulsive interparticle potential, plus 

£SJ ' gravity. This leads to a minimum lengthscale for bound objects, and to superfluidity. 

Galactic dynamics may differ significantly from that of more generic dark matter in 
J> ■ not unwelcome ways, especially in the core. Such particles can be realized as quanta 

of a relativistic massive scalar field with a quartic self-interaction. At high densities, 
, the equation of state has the same form as that of an ideal relativistic gas despite the 

interactions. If the nonrelativistic lengthscale is of order a kiloparsec, then the energy 
\ density in these particles was comparable to that of photons at early times, but small 

enough to avoid conflict with primordial nucleosynthesis. 

^ ■ 
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?— ( ' 1. Introduction 

^— > ■ 
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At the time of writing, the combination ("ACDM ") of cold dark matter, a Harrison-Zeldovich 
. !_< \ spectrum of initial fluctuations, and a small cosmological constant appears to be in good agreement 

with many observational constraints: the angular power spectrum of the temperature of the cosmic 
microwave background, the large-scale distribution of galaxies and clusters, and the growth of 
present-day structures by gravitational instability (Bahcall et al. 1999, and references therein). 
Recently, however, it has become clear that ACDM overpredicts structure on small scales. N-body 
and hydrodynamic simulations of galaxy formation within the ACDM model predict that the dark 
halos of galaxies should have singular cores (Navarro, Frenk, & White 1996; Navarro & Steinmetz 
2000), contrary to observation (Carignan 1985; Swaters, Madore, & Trewhella 1999). A separate 
though related problem is that late accretion of dark matter clumps or satellite galaxies is likely to 
shred the disks of spirals (Toth & Ostriker 1992; Klypin et al. 1999). 

Clustering on small scales could be suppressed by an upper limit to the phase-space density of 
the dark-matter particles due either to thermal entropy or, if they are fermions, degeneracy pressure; 
the particle mass would have to lie in the range 10 2 eV < m < 10 3 eV (Hogan & Dalcanton 2000). 

Another proposal is that the dark-matter particles are heavy (m 3> GeV) but self-interacting 
(Spergel & Steinhardt 1999). Since the simulations tend to produce a dark-matter velocity disper- 



- 2 - 



sion that decreases towards the center of the cusp, occasional collisions "heat" the central regions 
and produce a finite-density core, at least temporarily. Once the inner regions become roughly 
isothermal, the core will recollapse; the collision cross section must be chosen so that this does not 
happen within the age of the galaxy This proposal seems unlikely to alleviate the second problem, 
viz. the damage to galactic disks. 

This paper proposes a third solution to the problem of small-scale power: a Bose-Einstein 
condensate of dark matter particles, hence similar to the axion, but interacting via a repulsive 
potential of finite range. I show in §2 that cores would have a minimum size independent of their 
density or mass, and that the dark matter would behave as a superfluid. Such particles arise fairly 
naturally as quanta of a self- interacting scalar field (§3). For plausible choices of the minimum core 
size, the interaction energy per particle would have been comparable to the rest mass somewhat 
before the universe became matter-dominated, and at earlier times it would have a relativistic 
equation of state (p « p/3), hence slightly increasing the effective number of degrees of freedom in 
the radiation field, but not enough to violate the constraints of primordial nucleosynthesis. With 
respect to the standard cold-dark matter spectrum, density fluctuations in the linear regime would 
be suppressed on comoving scales less than a few megaparsecs (§4). 

While this work was being written up, I became aware that my colleague P. J. E. Peebles 
has been working along similar lines (Peebles & Vilenkin 1999; Peebles 1999, 2000). There is also 
closely related earlier work (Tkachev 1985, 1991). This seems to be a genuine case of convergent 
evolution, and the fact that independent groups arrived at similar results is perhaps reassuring in 
such a wildly speculative domain. I started with a nonrelativistic many-body view of these particles 
(§2) and then sought a relativistic framework for them (§3), whereas Peebles & Vilenkin (1999) 
seem to have proceeded in the opposite direction. With further development, the nonrelativistic 
but fully quantum-mechanical viewpoint may prove useful in studying the galactic dynamics of this 
form of dark matter. 



2. Minimum dark-matter core radius 

Suppose that nonrelativistic bosons interact via a a two-particle potential U(ri — r-2) of finite 
range. The potential energy of N such bosons in a common single-particle momentum state i/>(r) = 
V>-i/2gip-r in volume V is 

where U(p') is the fourier transform of U(r'), whose range is assumed to be small compared to the 
linear dimensions of V. If N, V — > 00 at fixed number density n = N/V, then the potential energy 
per unit volume is w(n) = n 2 U (0)/2. For the moment, w(n)/n <C mc 2 , the boson rest-mass energy. 

Macroscopically, one has a polytropic gas of adiabatic index 7 = 2; in other words, the pressure 
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is related to the mass density p ps ran by 

p = Kp 2 , K = U(0)/2m 2 . (2) 

If the gas is self-gravitating, there is a well-known spherical equilibrium with the density profile 
(Chandrasekhar 1939): 

. . . . sin(r/a) / K , . 

pW = K o,_LAi, „ = v /_. (3) 

The radius of the sphere, 7ra, is independent of the central density p(0), which determines the total 
mass, 47r 2 a 3 p(0). Dark-matter halos do not have the profile (3), but this does not rule out the 
model. If the particles are not all in the same momentum state, then their relative motions make 
an additional contribution to the pressure, which allows the halo to have a power-law density profile 
outside the core. Indeed, axionic dark matter is usually assumed to be a Bose-Einstein condensate 
like the one considered here but without the repulsive interaction, so the pressure support of dark 
halos in that model is due entirely to relative motions. Still, since (2) gives the minimum pressure, 
(3) is the most compact possible equilibrium. The conventional definition of the the core radius in 
terms of the central density and its second derivative is r c = [— 3p(0) / p" (O)] 1 ^ 2 , thus r c min = 3a. 

The nonrelativistic approximation breaks down when the pressure (2) is comparable to the 
rest mass density. This happens when p rs c 2 /Ga 2 = p re \. If this bosonic dark matter dominates 
the total present-day mass density, then the mean mass density was equal to p IC \ at redshift 

1 + z <« - ( w) V3 " 2 - lx 105 n ^*^ r ^' (4) 

where £1 = Sl/0.3, h^o = Ho/ (50 km s _1 Mpc -1 ), and r c k pc = fc,min/kpc. The question how 
the dark matter behaves at higher redshifts is deferred to §§3-4. It suffices for now that 1 + z re \ 
is comfortably larger than the redshift of matter-radiation equality (Peebles 1993): 1 + z eq = 
1.8 x 10 3 flo. 3 h 2 . 

The interaction makes the gas a superfluid. If one particle is removed from the Bose-Einstein 
condensate and put into a single-particle state with momentum p + q 7^ p, then the potential 
energy (1) is replaced by 



N ~ 1 

W N -t + — — 



j7(o) + r;(q) =w N + ^^u(d). 



The first term on the left is the interaction of the N — 1 particles in the condensate with one another, 
and the second is the interaction of the condensate with the extracted particle; the piece involving 
U (q) is the exchange energy resulting from symmetrization of the iV-particle wavefunction. I 
assume that the range of U(r) is sufficiently short so that U(q) ~ U(0). The energetic penalty for 
removing a particle from the condensate is then approximately equal to the potential energy per 
particle pair, nU(0). Thus if the condensate streams past an obstacle (an external potential) at 
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speed v, scattering out of the condensate is impossible if the kinetic energy per particle is less than 
this energy penalty, i.e. if the relative velocity 

v < ^2nU{0)/m = v crit {n). (5) 

Similarly, when two condensates of density m 2 and momenta pi t 2 stream through one another, 
dissipation occurs only if |pi — P2I > V2mv cr i t (ni + 712). This is not to say that the two streams do 
not interact, but rather that they interact only through the mean-field energy niri2U(0) per unit 
volume. Inside a core supported mainly by the repulsive interaction [eq. (3)], v CT n is comparable to 
the virial velocity. 

It will be interesting to study whether superfluid dark matter would have any distinctive con- 
sequences for galactic dynamics other than the minimum core size. Attention naturally focuses 
on dissipative processes, such as dynamical friction: i.e., irreversible transfer of energy and mo- 
mentum between the dark and baryonic matter via by their gravitational interaction (cf. Binney 
& Tremaine 1987). In collisionless systems, dynamical friction involves upon single-particle reso- 
nances (e.g. Tremaine & Weinberg 1984), much like Landau damping in plasmas. As long as the 
condition (5) is satisfied, however, a perturbing gravitational potential interacts coherently with 
the condensate, and all particles have the same resonant frequencies because they share a common 
macroscopic wavefunction. Thus for example, a rotating galactic bar may experience little drag 
against the dark matter; this may circumvent an important argument against dense dark halos 
in barred spirals (Debattista & Sellwood 1998). The question will require a quantitative analysis, 
however, because even in the innermost parts parts of the galaxy, not all of the dark matter will 
be in the condensate. 



3. Relativistic era 

Dark matter with the properties described in §2 arise as quanta of a self-interacting relativistic 
scalar field <p with lagrangian density 

c = -V=g (^g^d^dA + v(</>)) . (6) 

Without loss of generality, the minimum of V(4>) occurs at (p = 0, and V(4>) = m 2 (j) 2 /2 + 
(higher powers). Potentials of the form 

y(0) = imV + ^ 4 (7) 

are of particular interest, though one might want to add a constant V(0) = K/SttG to produce a 
present-day cosmological constant. In lowest-order perturbation theory, the interaction energy of 
a state |^tv(0)) consisting of N quanta at rest in volume V is, in Minkowski space, 

J d\ (<M0)| : ^ 4 (r, t) : |^(0)) = ^ ^= 1} , (8) 
V 
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and therefore U(0) = 3K/m 2 . 

Semiclassical methods give the same result, which is important because they are not restricted 
to perturbation theory. Thus if <ft were a spatially uniform classical field, then C could be regarded 
as the lagrangian of a one-dimensional oscillator with an explicit time dependence via the metric 
9iiv diag(— 1, a 2 , a 2 , a 2 ) in an Einstein-de Sitter universe. The momentum conjugate to (ft is 
w = a 3 (ft, the hamiltonian is 

2 

H = ^ + a 3 V(<ft), 

and the action in the oscillator is given by an integral over one complete cycle: 



If,, a 



3 



1 = — i> wdcft = — = <b y/a~ 3 H - V{(ft) d<ft. (9) 
2vr J vrV2 J 

Semiclassically, X becomes the number of quanta per comoving volume, na 3 , while Ji becomes the 
energy per comoving volume, pa 3 . (In this section, p will be the total energy density, not the 
rest-mass density ran of §2.) The definition (9) makes sense only when the oscillation frequency 
to = {&H/dI) a is much larger than the current Hubble expansion rate a/a, in which case X is an 
adiabatic invariant and hence na 3 is conserved. By direct expansion of the quadrature (9) to first 
order in k and inversion of series, one has 

3k 

hence p = mn -\ ^ n 2 + 0(K 2 n 3 ), (10) 

2m z 

in agreement with previous results for the nonrelativistic (small-n) regime. In the opposite limit 
of large n, the quadrature (9) is dominated by (ft S> tti/^/k; neglecting the mass term in V((ft), one 
has 

p « 3 4 / Vr- 8 / 3 (l/4) «V3 n 4/3 „ 1 377 K i/3 n 4/ 3; (n) 

as if this were a noninteracting relativistic gas: p = —d(pV)/dV = p/3. Eq. (9) can be evaluated 
to an exact expression for n(p) in terms of complete elliptic integrals. 

We are now in a position to estimate the mass m and average number density n{z) of these 
quanta. From eqs. (3) & (10), it follows that the minimum core radius r C)in i n = 3a depends only 
upon m A /K and fundamental constants, so 



- 2 =(^S^V /4 -10.7^^eV. (12) 



c.min 



Apart from the dimensionless coupling k, this is the geometric mean of the Planck mass and the 
mass whose Compton wavelength is 27rr Cim j n . Furthermore, if this form of dark matter dominates 
the mass density today, then 

n{z) = 5^(1 + z) « 74. ^'/^fio,^ (1 + *f cm- 3 . (13) 
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Prior to the redshift (4) when particles followed the relativistic equation of state (11), they would 
have contributed a constant fraction of the total energy density, equivalent to an increase 

AiV^0.14r c 2/ k 3 pc (Oo.3^o) 4/3 (14) 

in the number of effectively massless neutrinos (assuming N v ~ 3), which is compatible with the 
constraints from primordial nucleosynthesis (Olive & Thomas 1999). 

4. DISCUSSION 

We have seen that small-scale structure can be suppressed even if the dark matter is completely 
cold and bosonic, provided that it has a repulsive interaction. At first blush, the idea seems less 
natural than the alternatives — warm or degenerate fermionic dark matter — which have been much 
more widely discussed. I are not aware of a strong particle-physics motivation for matter with these 
properties. 

Nevertheless, in working through the consequences of the basic idea, one is intrigued by some 
satisfying coincidences. 

(i) From a nonrelativistic viewpoint (§2), the equation of state (2) results from a generic two- 
body interaction of finite range among massive bosons; relativistically, it emerges from the 
simplest nonlinear field theory (6)-(7). 

(ii) The nonrelativistic equation of state implies a characteristic lengthscale and a minimum 
dark-halo core radius. If this lengthscale is of order a kiloparsec, as the observations suggest 
(Hogan Sz Dalcanton 2000), then the dark matter began to be nonrelativistic at the lowest 
possible redshift that growth of structure would permit, viz. z rc \ ~ z ec[ . The result (14) that 
the energy density in these hypothetical quanta would have been comparable to the energy 
density in photons at early times is really the same coincidence. Both are independent of n 
and m, because the relevant combination of these quantities is already fiexed by r Cim i n & 17. 

With regard to the second point, Peebles (2000) has estimated that the model may be a little 
too successful at suppressing small-scale power during the linear regime. Density fluctuations that 
come within the horizon before z eq not only do not grow, but actually decay, until their physical 
size is larger than the "Jeans length" a. He finds that this constraint is marginally inconsistent 
with the quartic model unless r Cjin i n < 0.5 kpc. Pending more precise observations, however, one 
may be impressed that the model marginally survives this test without appeal to an adjustable 
parameter. 1 



1 Peebles notes that slightly sub-quartic potentials V((f>) = m 2 <j) 2 /2 + K,\4>\ q with q « 3.7 would fit this constraint 
more comfortably. But one would have to sacrifice (i). 
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